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ABSTRACT 

The accretion induced neutron star magnetic field evolution is studied through considering the 
accretion flow to drag the field lines aside and dilute the polar field strength, and as a result the 
equatorial field strength increases, which is buried inside the crust on account of the accretion 
induced global compression of star crust. The main conclusions of model are as follows: (i) the 
polar field decays with increasing the accreted mass; (ii) The bottom magnetic field strength 
of about 10® G can occur when neutron star magnetosphere radius approaches the star radius, 
and it depends on the accretion rate as Af (iii) The neutron star magnetosphere radius 
decreases with accretion until it reaches the star radius, and its evolution is little influenced 
by the initial field and the accretion rate after accreting ^ O.OIA/©, which implies that the 
magnetosphere radii of neutron stars in LMXBs would be homogeneous if they accreted the 
comparable masses. As an extension, the physical effects of the possible strong magnetic zone 
in the X-ray neutron stars and recycled pulsars are discussed. Moreover, the strong magnetic 
fields in the binary pulsars PSR 1831-00 and PSR 1718-19 after accreting about half solar 
mass in the binary accretion phase, 8.7 x 10^" G and 1.28 x 10^^ G, respectively, can be 
explained through considering the incomplete frozen flow in the polar zone. As a model's 
expectation, the existence of the low magnetic field 3 x 10^ G) neutron stars or millisecond 
pulsars is suggested. 
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1 INTRODUCTION 

Recently, the newly discovered double pulsars J0737-3039A and 
J0737-3039B in the binary system (Burgay et al 2004; Lyne et al 
2004; van den Heuvel 2004; Lorimer 2004) have shown that the 
millisecond pulsar (MSP) possesses the low field and short spin 
period (B = 7x 10^ G, P=22.7 milliseconds) and the normal pulsar 
possesses high field and long spin period (B = Qx 10^^ G, P=2.77 
seconds), which convinces our previous point of view that MSP is 
recycled in the binary system where the accreted matter weakened 
its magnetic field and accelerated its spin (see, e.g., van den Heuvel 
2004). Moreover, with the launch of RXTE satellite, the accretion- 
powered X-ray pulsar SAX J 1808.4-3658 (P=2.49 milliseconds), 
as well as other five examples (see, e.g., Chakrabarty 2004; van der 
Klis 2004; Wijnands & van der Klis 1998), and 11 type-I X-ray 
burst frequency, confirmed as stellar spin frequency, sources have 
been discovered (see, e.g., Muno 2004; Chakrabarty 2004; van der 
Klis 2004, 2000), which exhibits the strong evidence that neutron 
star (NS) in low mass X-ray binary (LMXB) is the progenitor of 
MSP with spin frequency ~ 400 Hz and magnetic field strength 



~ 10* ^ G because of the accretion (see, e.g., Wijnands & van der 
Klis 1998; Chakrabarty 2004). 

The origin, structure and evolution of NS magnetic field is still 
an open problem (see, e.g., Bhattacharya & van den Heuvel 1991; 
Phinney & Kulkarni 1994; Bhattacharya & Srinivasen 1995; Colpi 
et al 2001; Melatos & Phinney 2001; Gumming 2004; Blandford et 
al 1983; Blondin & Freese 1986). In the early stage of discovery 
of NS, it is believed that NS magnetic field decays due to Ohmic 
dissipation in the crust. But later calculations of Ohmic dissipation 
suggest that isolated NS magnetic field may not decay significantly 
(Sang & Chanmugam 1987), if the field occupies the entire crust, 
and remains large for more than Hubble time. Although many in- 
structive mechanisms on the isolated NS magnetic decay have been 
proposed, such as, the crustal plate tectonics model by Ruderman 
(1991) and Chen & Ruderman 1993, the Hall-drift dominated field 
evolution (see, e.g., Rheinhardt et al 2004; Jones 2004; Hollerbach 
& Ruediger 2002; Geppert & Rheinhardt 2002; Rheinhardt & Gep- 
pert 2002; Naito & Kojima 1994), the spin-evolution induced mag- 
netic field decay (see, e.g., Konar & Bhattacharya 1997, 1999; Ding 
et al 1993; Jahan Miri & Bhattacharya 1994; Ruderman et al 1998), 
etc., there has not yet been a commonly accepted idea on such is- 
sue. 
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paid much attention until 1980's when Taam and van den Heuvel 
(1986) presented the NS magnetic field evolution associated with 
the accretion in the binary system, based on which Shibazaki et 
al (1989) concluded the simple empirical formula of the field de- 
duction versus the accreted mass. However, it is indicated that the 
significant decay of magnetic field is achieved only if the neutron 
star experiences the interacting binary (see, e.g., Verbunt & van den 
Heuvel 1995; van den Heuvel 1995 and references therein). Fur- 
thermore, van den Heuvel and Bitzaraki (1995), from the statistical 
analysis of 24 binary radio pulsars with nearly circular orbits and 
low mass companions, discovered a clear correlation between spin 
period and orbital period, as well as the magnetic field and orbital 
period. These relations strongly suggest that the increased amount 
of accreted mass leads to the decay of NS polar magnetic field, and 
the 'bottom' field strength ~ 10* G is also implied, which means 
that NS magnetic field evolution will stop at some minimum value, 
~ 10** G (see also, e.g., Burderi et al 1996; Burderi & D'Amico 
1997). 

In order to understand the accretion induced field decay, some 
suggestions and models have been proposed. Based on the ac- 
cretion buried or screened NS magnetic field guess by Bisnovati- 
Kogan and Romberg (1974), Zhang et al (1994, 2000) proposed 
the ferromagnetic crust screen model to interpret the simple inverse 
correlation between the field deduction and the accreted mass, as 
declaimed by Shibazaki et al (1989). Moreover, the accretion in- 
duced flow and thermal effects to account for the speed-up Ohmic 
dissipation of NS crust currents are also proposed and studied by 
several authors (see, e.g., Romani 1990; Geppert & Urpin 1994; 
Urpin & Geppert 1995). Later on, the diamagnetic screening of NS 
magnetic field by freshly accreted material is suggested (see e.g. 
Lovelace et al. 2005), which would be effective above a critical ac- 
cretion rate Al of a few percent of the Eddington rate (see, e.g.. 
Gumming et al 2001). More recently, it is calculated that the polar 
cap widening or the equator-ward hydro-magnetic spreading by the 
accreted material flowing accounts for the NS magnetic field de- 
cay (see, e.g., Konar & Choudhury 2004, references therein; Payne 
2005; Payne & Melatos 2004; Melatos & Phinney 2001). Further- 
more, the magnetic burial mechanism has been studied in detail an- 
alytically and numerically by Payne & Melatos (2004) and Payne 
(2005) under the assumptions of no instabilities or Ohmic diffusion, 
where the structure of highly distorted magnetic field and hence the 
magnetic moment as a function of accreted mass is computed. They 
find that the deduction of the magnetic moment is inversely related 
to how much mass accreted and about ~ W'^Mqis required to 
reduce significantly the magnetic dipole moment. 

In this paper, we follow the idea, firstly proposed by van den 
Heuvel and Bitzaraki (1995), and also by Burderi and D'Amico 
(1997), that the accretion matter is channeled by the strong mag- 
netic field initially into the two polar patches, corresponding to NS 
magnetosphere radius of several hundreds of stellar radius, pushes 
the field lines aside and thus dilutes the polar field strength due to 
the flow of accreted materials from the polar to the equator. The 
bottom field should be reached when the accretion becomes all 
over star, which corresponds to that the NS magnetosphere radius 
matches the star radius, and in turn it infers the bottom magnetic 
field to be ~ 10* G. 

The paper is organized as follows: in the next section the 
model is described and the magnetic field evolution equation is de- 
rived. Section 3 outlines the applications of the model, including 
how the NS magnetic field and magnetosphere evolve with the ac- 



discussions are drawn in Section 4, where we summarize the con- 
clusions of model and the observational consequences. 



2 DESCRIPTION AND CALCULATION OF THE MODEL 
2.1 Dilution of the polar magnetic flux 

For the spherical magnetic slab geometrical structure of accreted 
NS illustrated in Fig. 1, we assume that the magnetic field lines 
are anchored, completely frozen or incompletely frozen, in the en- 
tire NS crust with constant average mass density. The introduc- 
tion of the spherical topological magnetic polar cap geometry is 
on account of its important application in describing the kHz QPO 
mechanism (Zhang 2004), where the spherical magnetic polar cap 
geometry in Eq.J^ is used and the approximated plane magnetic 
slab (Cheng & Zhang 1998) is inadequate when the magnetosphere 
is close to the star surface. Moreover, in the work by Cheng & 
Zhang (1998) the rate of change of the magnetic polar cap area is 
roughly performed (see Eq.(5) of that paper), and it should be de- 
rived from the magnetic flux conservation as described below. Ini- 
tially, the magnetic field is sufficient strong, and the accreted matter 
will be channeled into the polar patches by the field lines. The cu- 
mulated accreted matter will be compressed into the crust region 
through the polar patch or through the surface flow on account of 
the plasma instabilities (see, e.g., Litwin et al. 2001). However, as a 
preliminary exploration, we neglect these instability effects on the 
polar field deduction, and concentrate on the description of the field 
evolution by the assumed frozen-in crust MHD motions. We sup- 
pose the compressed accreted matter to arise the expansion of the 
magnetic polar zone into two directions, downward and equator- 
ward. Under the condition of incompressible fluid approximation 
and the constant crustal volume assumption, as well as the aver- 
aged homogenous NS mass density p ~ W^'^g/cm^, where the 
Ohmic dissipation time scale of field is longer than the flow time 
scale for sub-Eddington accretion rates (see, e.g., Brown & Bild- 
sten 1998; Geppert et al 1999), the magnetic line frozen motion is 
assumed and the incomplete frozen motion may appear if the un- 
known instabilities are taken into effects (see, e.g., Melrose 1986). 
Therefore, in a particular duration St of accretion with the mass ac- 
creting rate A/, the piled accreted mass SM = MSt will arise the 
expansion of the volume of the magnetic polar zone in the crust, by 
virtue of the mass flux conservation, 

{SM/2 - pApSH)^ = pHSAp , (I) 

where Ap is the surface area of the magnetic polar zone defined in 
Eq.j4j- The introduced parameter ^ (0 ^ ^ ^ 1) is an efficiency 
parameter to express the incomplete magnetic line frozen flow by 
the plasma instability. In two extremely cases, ^ = 1 represents 
that the completely frozen field lines totally drift with the moving 
mass elements, and ^ = represents that the completely leaking 
of plasma makes the field lines not drift with the moving mass ele- 
ments. The physical meaning of Eq.Q is clear, ^ = 1 for instance, 
that the r.h.s of Eq.Q is the mass flux toward the equatorial direc- 
tion and the l.h.s of Eq.Q is the mass accreted in one polar cap 
(SM/2) minus the mass sinked into the NS core. H is the thickness 
of the crust and SHfH is the fraction of the thickness dissolving 
into the core, defined by 

SH SM MSt 
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where Mcr = 4:TtR^pH is the crust mass. For the standard neu- 
tron star structure model, the NS crust mass is only a small por- 
tion of the entire NS mass, which is also dependent of the equation 
of state (EOS) of matter inside NS (see, e.g., Lattimer & Prakash 
2004; Baym & Pethick 1975), so M^r/M ~ 3H/R ~ 10% is 
approximately implied. The expansion of the polar zone will di- 
lute the magnetic flux density if the magnetic flux conservation is 
preserved, that is. 



5{BAp) = 



SAp = 



A ^ 



(3) 



and the area Ap of the accretion polar patch can be accurately 
expressed by the angle 9, which is subtended at the stellar surface 
between the last closed field line of the NS magnetosphere and the 
polar axis. 



Ap — 



2tvR^ sin e'de' = 27ri?^ ( 1 - cos f 



(4) 



The polar field line can be described by the polar angle 9r and the 
position r that satisfies (sin'^ 9r) /r=constant (see, e.g. Shapiro & 
Teukolsky 1983, p453), therefore, at the stellar surface, in terms of 
the NS magnetosphere radius Rm, 



R 

Rm 



(5) 



If 6 is small as a case of large magnetosphere radius, we have the 
approximated relation Ap = tvR^ sin^ 9 = ttR''/Rm, as conven- 
tionally used (see, e.g. Shapiro & Teukolsky 1983, p453). 
Rm is related to the Alfven radius Ra by 



Rm = <j)RA , 



Ra = 3.2 X 10'*(cm)Mj"/'^V30^"i~^'^'^ 



(6) 



(7) 



where (pis a parameter of about 0.5 (see, e.g., Ghosh & Lamb 1979; 
Shapiro & Teukolsky 1983), but is model dependent (see, e.g., Li 
& Wang 1999). m = M/Mq is the NS mass M in units of solar 
mass, and Alir and /iao are the accretion rate in units of 10^^ g/s 
and the magnetic moment in units of 10^° Gcm^, respectively. 

2.2 The buried magnetic flux in the equatorial zone 

With the proceeding of the accretion, the magnetic flux in the polar 
zone will be diluted and the magnetic flux in the equatorial zone 
would be increased because of the global conservation of the mag- 
netic flux. However, the new constructed magnetosphere radius is 
still determined by the new constructed polar magnetic field be- 
cause the increased equatorial magnetic flux is compressed into the 
star crust by the global radial accretion flow. Therefore the equato- 
rial strong magnetic field lines are curved and buried inside crust 
and might be not seen during the accretion phase with the suffi- 
cient accretion rate, but it may be seen in the post-accretion phase 
on account of the Ohmic diffusion. We will try to confirm the above 
claim through comparing the magnetic flux drifted from the polar 
zone to the equatorial zone and the magnetic flux compressed into 
the crust in the equatorial zone. 

We substitute Eq.j2j into Eq.Q, and then obtain 



(27ri?^ - Aj_ 



, ^SH SAp 



HSt 



St 



(8) 



we can obtain the ratio between the field line drift velocities in the 
latitude direction and in the radial direction. 



^Rcos9 
Hsm9 



(9) 



with the radial velocity = SH/St — M /AnR^p and the latitude 
velocity ve = RS9/St. 

The polar magnetic flux variation rate expelled from the polar 
zone with the approximated radius R sin 9 in the latitude direction 
is 

^2nRBsm9ve , (10) 

with $i„ — '^'^j" , and the equatorial magnetic flux rate com- 
pressed into the crust in the radial direction is 

^out~2nRBvro , (11) 

^TT^ and Vio = M /^ttR^ po, where po is a charac- 



with $o„t = 

teristic mass density at where the matter pressure controls the mag- 
netic pressure. If the non-relativistic neutron in the crust is taken 
into account, the matter pressure is expressed to be Pm = kp''^^'^ 
with k=5.38 x 10'-' in c.g.s. units (Shapiro & Teukolsky 1983). And 
the magnetic pressure Pb is calculated by the possible maximum 
permitted magnetic field strength Beqmax ~ 10^^ Gauss during 
the accretion compression, exceeding over which the magnetic en- 
ergy will be released by the soft gamma ray burst triggered by the 
super-strong magnetic induced crust cracking (see, e.g., Thompson 
& Duncan 1995), therefore we have. 



Bi 



Stt = fcpo 



5/3 



or 



Po ~ 10"(<?/cm')(Be,™a./10'')'/' 



(12) 



(13) 



The ratio between the flux variations in the polar zone and in the 
equatorial zone can be given if the condition of Eq.j9} is consid- 
ered. 



_ V0sai9 _ £^Rpo cos 9 

^out «ro Hp 



(14) 



0.0l^cos9{B,q^ajW'^''f^^ < 1 



and we find that this ratio is less than unity for the possible param- 
eters, which means that the magnetic flux expelled into the equato- 
rial zone through the latitude direction is much less than that com- 
pressed into the crust through the radial direction. Therefore, we 
can conclude that the magnetic field strength at the magnetosphere- 
disk boundary will be dominated by the field lines from the polar 
cap zone, or the magnetosphere radius Rm during accretion will be 
determined by the polar field. The effective magnetic moment of 
star decreases because of the motions of field lines, dragged firstly 
to the equator zone from the polar zone through the latitude mo- 
tion and then compressed into the crust through the radial motion. 
Physically, the open field line in the polar zone will be influenced 
by the latitude flow and little influenced by the radial flow, however 
the closed field line in the equatorial zone will be dominated by the 
radial flow. If there exists only the radial flow and no latitude flow, 
as a case of the magnetosphere reaching the star surface, the po- 
lar field will be little changed because the net motion of field lines 
in the latitude direction disappear Therefore, with the proceeding 
of the accretion, the polar field decays and the newly constructed 
magnetosphere radius is still determined by the decayed polar field 
strength, or in mathematical terminology the field in the Eqs.JU 
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Furthermore, we stress that the above arguments for the field 
hne motion are roughly vahd for the purpose of the phenomenolog- 
ical model, however the detail descriptions of the accretion flow, 
together with the magnetic induction equations, are needed in pre- 
senting the NS magnetic structure, which will be the subsequent 
work elsewhere. 



2.3 The polar magnetic field evolution 

Furthermore, substituting Eqs.j2}, J3}™d ^4) into Eq.Q, we obtain 
the magnetic field evolutionary equation as follows. 



ApSB 



(27ri?2 _ Ap)B 



or, 



V-Rm - R)SB _ -iMSt 



(15) 



(16) 



Solving Eq. <15t or Eq. <16> with the initial condition B(t=0)=Bo 
the analytic magnetic field evolution solution is obtained. 



B 



Bi 



{l-[C/exp(y)-l]2}V4 



(17) 



with , = AM 

2 _ (|,)4/r 



Mt, C ^ 1 + \J\-xi ~ 2 and 

Xo~ = (^j"' ' = -R/Rmo, where Rmo is the initial magneto- 
sphere radius and -B/ is the bottom magnetic field defined by the 
NS magnetosphere radius matching the star radius, i.e., Rm{B{) = 
R, therefore, 



Bf = 1.32 X 10* (G) ( 



M 



(18) 



where M^d is the Eddington accretion rate and Rq is the NS ra- 
dius in units of 10® cm. For the reason of simplicity, the usual NS 
parameter values are set in the later calculations and discussions, 
such as, m=1.4, i?6 = 1-5, cj) = 0.5 and = O.2M0. 



3 THE APPLICATIONS OF THE MODEL 

3.1 Inverse correlation between the magnetic field and the 
accreted mass 

In the complete magnetic frozen case ^ = 1 (for the incomplete 
frozen case, we just take the replacement Mcr/C ~^ Mcr), the 
magnetic field evolution with the accretion in Ea. <17> can be ap- 
proximately simplified in the following form by the Taylor serial 
expansion with the condition y ^ 1 or AM ^ 3.5Mcr, 

Bo 



B = 



Q I 4AM \ 7/4 
^ 7m D ' 



(19) 



where niB ~ Q.^Xo Mcr = 0.5(-R/RMo)Mcr. For the usual value 
of bottom field Bf ~ 10* G and initial field Bo ~ 10^^ G, we 
obtain ms ~ Mq. For the high mass X-ray binary (HMXB) 
or the low mass X-ray binary (LMXB) in the early stage, a little 
mass is accreted, AA/ ~ 10~"'Mq <C ms, then Ea. <19t gives, 

Bo 



B 



1 -(_ 



(20) 



fitting the observational data given by Taam and van den Heuvel 
(1986). Shibazaki et al (1989) find that the best fitted mass constant 
is TOfl = 10~® - 10~^ M0when comparing the theoretical curve of 
the magnetic field versus the spin period to the observational data. 
Moreover, the similar mass value of significantly reducing the mag- 
netic dipole moment is also obtained by Payne &. Melatos (2004) 
and Payne (2005) in their self-consistent analytic and numerical 
calculations of the accretion buried NS polar magnetic field. 

If ms <C AM <^ 3.5Mcr, then we have the following ap- 
proximation from Eq. <17> . 



(21) 



Therefore, Eq. <21> implies that the influence of the initial magnetic 
field on the magnetic evolution disappears at this stage, and the NS 
magnetic field is scaled by the bottom field Bi. Furthermore, with 
increasing of the accreted mass AM Mcr, the NS magnetic evo- 
lution will enter into the "bottom state", and B —> B{ is obtained. 
As a detail illustration, the evolution of NS magnetic field versus 
the accreted mass with various parameter conditions is plotted in 
Fig. 2, where we find that the influence by the initial field on the 
field evolution exists when AM <C O.OIM© and disappears when 
AM > O.OIMq. Finally, the magnetic field evolution will enter 
into the "bottom state" when AM Mcr ~ 0.2Mq. 



3.2 The bottom magnetic field and its correlation to the 
X-ray luminosity 

The conception of bottom field, the possibly arrived minimum 
field of accreted NS, is evidently proposed by van den Heuvel and 
Bitzaraki (1995) from the analysis of the magnetic fields of the mil- 
lisecond pulsars in the binary systems versus the estimated accreted 
masses, which is explained that the accreted matter pushes and di- 
lutes the polar field lines if the accretion in the channel way and 
there is no net flow drag effect on the field lines if the spherical 
accretion in the random way all over star begins. From the point of 
view of the accretion geometry, the bottom field is reached when 
the accretion becomes isotropic all over star or the channeled ac- 
cretion disappears completely (see also, e.g., Burderi et al 1996; 
Burderi & D' Amico 1997). 

Mathematically, as a minimum field, the bottom magnetic 
field can also be determined by the condition of vanishing the field 
variation respect to the accreted mass in Eq. <17> . i.e., 



SB 
5{AM) 



= 0, 



(22) 



which gives the magnetic field minimum value to be the bottom 
field B{ when the following condition is satisfied. 



C/exp(j/)-l = 0, 



(23) 



which is just the same form as the empirical formula of the accre- 



which can give a critical accreted mass, namely, AM — 
3.5Log{C)Mcr ~ O.5(M0)(Mcr/O.2), after which the bottom 
field will be preserved whether or not how much extra mass ac- 
creted. From the statistics of 24 binary pulsars by van den Heuvel 
and Bitzaraki (1995), this critical mass is approximately a fraction 
of ~ O.7M0. 

In Fig.2, the influence of the accretion rate (luminosity) on the 
magnetic field is clear that the Eddington luminosity Lag (Lse = 
O.OlLss) corresponds to the bottom field 3 x 10* G (3 x 10^ G). 

The relation between the X-ray luminosity and the mas- 
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Hasinger and van der Klis (1989), and they found that the Z (Atoll) 
sources with the Eddington luminosity (~ 1% Eddington luminos- 
ity) possess the strong (weak) magnetic fields. However, we stress 
here that the mass accretion rate cannot be exactly inferred from 
the observations, therefore we conclude that the individual Atoll 
source possesses a weaker field than that of Z source only when 
both sources accretes the similar masses and the mass accretion 
rate of Atoll source is lower than that of Z source. 

Moreover, the theoretical analysis also implies the propor- 
tional correlation between the magnetic field and the X-ray lumi- 
nosity (see, e.g., Psaltis & Lamb 1998; Compana 2000). However, 

1/2 

the magnetic field and luminosity relation of NS/LMXB B oc LJ 
was also proposed by White and Zhang (1997) when analyzing the 
kHz QPO data of the ten LMXB samples under the assumption 
that the fast spinning neutron stars are close to the spin equilib- 
rium state. In reality, for the specific magnetosphere radius Rm 
of LMXB, if scaled as Rm oc (B^/M)^'"' ~ R, the bot- 
tom field should be proportionally related to the accretion rate as 
Bf oc Afi/^ oc lV^. 

3.3 The effects of parameter f 

If ^ <C 1, the flow dragging effects on the field lines will be low 
efficiency and the field decay will be slow. In other words, the small 
^ causes decreasing the effective accreted mass contributed to the 
field decay from Ea. <17> . The ^ parameter influence can be applica- 
ble to the binary pulsars PSR 1831-00 and PSR 1718-19, believed 
to possess strong magnetic field 8.7 x 10^° G and 1.28 x 10^^, 
respectively, after accreted about half solar mass in the binary ac- 
cretion phase (van den Heuvel & Bitzaraki 1995). Thus we may 
imagine that the progenitors of both pulsars experience the incom- 
plete frozen plasma flow in their magnetic polar caps, on account 
of the unknown instabilities, the mechanism of which is still un- 
clear. For these two systems, as drawn in Fig.2, if we set ^ ~ 0.03 
and ^ ~ 0.002, their polar magnetic field can decay to the present 
values from the assumed initial value ~ 5 x 10^^ G after accret- 
ing half solar mass. While, the introduction of ^ parameter helps 
us to interpret the sustained strong magnetic field NS after accret- 
ing ~ O.5M0 in systems PSR 1831-00 and PSR 1718-19 (van den 
Heuvel & Bitzaraki 1995), and we find a parameter other than the 
accreted mass to influence on the magnetic field decay, which was 
expected by Wijers (1997). Moreover, as another example, it was 
also pointed out by Verbunt et al (1990) that NS in 4U 1626-67 has 
accreted ~ O.IM0 without losing its magnetic field ~ lO^'^ G. 

3.4 The NS magnetosphere evolution 

The NS magnetosphere radius can be defined by the magnetic field 
through the condition Rm/R = (B/B/)*''^, 

Rm ^ 1 ^ exp(2y) 

R l-[C/exp(y)-l]2 2C[exp(y) - C/2] • 

If y <C; 1 or AM <^ 3.5Mcr, then we have the approximated cor- 
relation between the magnetosphere radius and the accreted mass, 

Rmo -Rm 1 ,,,, 

= °'' ^ R I 8AM • (25) 

^ rmg Rmo 7Mor 

Furthermore, if tub «C AM ^ 3.5Mcr, then we obtain 

AM , 



which implies that the magnetosphere radius evolution at this stage 
has nothing to do with any initial conditions and even the X-ray 
luminosity and only depends on the accreted mass. If the accreted 
masses are similar for the NSs in LMXBs, their magnetosphere 
radii would be homogeneous for both high luminosity Z sources 
and low luminosity Atoll sources. At the last stage of LMXB, if 
AM ^ Mcr, then we have Rm ^ R- The above analytical conclu- 
sions are illustrated in Fig. 3, where the magnetosphere evolution 
versus the accretion mass is plotted, and the accretion rate influence 
and the initial field influence are taken into account. We find that the 
magnetosphere radius is proportionally related to field strength and 
inversely related to the accretion rate at the initial evolution stage, 
however it has little correlation to the accretion rate (or luminosity) 
and the initial field when the the accretion mass exceeds over about 
O.OIA/q . In other words, in the late evolution stage, LMXB for in- 
stance, the NS magnetosphere radius is scaled by the NS radius, 
and only depends on the accreted mass. 

3.5 The strong magnetic field estimation in the equator zone 

With the accretion going, the magnetic field lines are dragged by 
the accretion flow from the polar cap toward the equator zone, at 
where they are buried inside the crust by the global compression 
of the accreted matter or decayed by the possible reconnection of 
the field lines in the crust. At the last evolution stage of LMXB, 
the global magnetic flux conservation will help us to estimate the 
buried equator magnetic field by assuming the initial magnetic flux 
to equal the final magnetic flux, 

2-kR'^Bo ~ 2TvR^Bi + 2nRHBeq . (27) 

For the usual condition Bf <^ Bo, the magnetic field in the equa- 
tor zone buried inside the crust can be estimated to be Beq — 
{R/H)Bo, which is much stronger than the initial field value. Well, 
it is remarked that this estimation of the equatorial magnetic field is 
based on the assumption of homogeneously distributed field lines 
inside the crust. On the detailed magnetic configuration geometry 
with the accretion, we refer to the analytical and numerical calcu- 
lations by Payne & Melatos (2004) and Payne (2005), where they 
compare the magnetic field line structures before/after accretion. 
As for the real NS mass density profile in the crust, the lower the 
mass density, the faster the accretion flow velocity. So the more 
field lines may be cumulated in the low mass density layer, where 
the field strength would be high but lower than B^qmax ~ 10^^ 
Gauss, because the fields as large as this would be unstable in the 
outer crust and might be able to break the crust at low enough den- 
sities (see, e.g. Thompson & Duncan 1995). The observational ev- 
idence to support the existence of the strong field may occur in the 
Type-I X-ray burst sources, which are due to unstable thermonu- 
clear burning of accreted hydrogen and helium on the NS surface 
to arise the modulation of X-rays at the star rotation, and it is be- 
lieved that the burst is firstly ignited in the hot spot and then spreads 
all over star (see, e.g., Strohmayer & Bildsten 2003; Gumming 
2004). Therefore, we associate the hot spot with the existence of 
strong magnetic area. Moreover, the accretion-powered X-ray pul- 
sar SAX J1808.4-3658 (Wijnands & van der Khs 1998; Wijnands 
et al 2003) hints the magnetically channelled accretion, which may 
also not exclude the existence of strong magnetic field or the ef- 
fects of multiple magnetic moment (see, e.g., Psaltis & Chakrabarty 
1999). As a plausible application, recently. Cutler (2002) points out 
that an internal strong magnetic field, while keeping the external 
dipole magnetic field low, may arise NS mass quadruple to produce 
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braking. As a result, NSs in various X-ray binaries and recycled 
millisecond pulsars could then be detectable by advanced gravita- 
tional wave interferometers. In addition, if the details between the 
accretion and polar field are taken into account, as pointed out by 
Melatos & Payne (2005), during accretion the NS magnetic field is 
compressed into a narrow belt at the magnetic equator by material 
spreading equator-ward from the polar cap. In turn, the compressed 
field confines the accreted material in a polar mountain, which is 
misaligned with the rotation axis in general, to produce the gravita- 
tional waves. At last, we stress that the magnetic field is not totally 
decayed but redistributed on account of the accretion flow, and the 
polar field decreases and the equator field increases but buried in- 
side the crust. The postaccretion diffusion of the suppressed field 
back to the NS surface may happen in the recycled pulsars as de- 
scribed by Young & Chanmugam (1995). Recent observations, both 
in the X-ray range of PSR 1821-14 (see, e.g., Becker et al. 2003), 
and in the radio range (see, e.g., Gil et al. 2002), seem to support the 
idea of the existence of strong small-scale magnetic field structures 
at the NS surface, and the further exploration of which has been 
paid attention recently by a number of authors (see, e.g., Geppert et 
al 2003; Uipin & Gil 2004). 



Payne & Melatos (2004) and Payne (2005), however the physical 
meaning of latter is at where the hydrostatic pressure at the base of 
accretion column overcomes the magnetic tension and the matter 
spreads over the stellar surface to drag the polar field lines toward 
the equator. In our model, tub represents the fraction of the crust 
mass scaled by the ratio of the stellar radius to the initial magneto- 
sphere radius, then the relation of both characteristic masses is still 
unclear. Finally, although our model has given the observational 
consistent results, it is a very simple and crude theoretical frame- 
work to account for the accretion induced polar field decay. We 
think that the basic physical idea is clear and has been proposed by 
a couple of people (see, e.g., van den Heuvel and Bitzaraki 1995; 
Burderi et al 1996; Burderi & D'Amico 1997) that the transfor- 
mation from the systematic channeled accretion to the spherical 
accretion yields the power to decrease the NS polar field strength. 
However, many physical details have been neglected in the model, 
such as, the detail of accretion flow and its reaction with the field 
lines, the plasma instability in the NS melted surface, the details of 
how the magnetic field is buried in the equatorial zone, the Ohmic 
dissipation, as well as the HaU-drift effects, etc., therefore the con- 
siderations of these will construct our future research explorations. 



4 DISCUSSIONS AND CONCLUSIONS 

In summary, the main conclusions of model are listed in the fol- 
lowing: (1) The NS polar magnetic field decays in the binary accre- 
tion phase, and its evolution experiences the following processes. 
Roughly speaking, B oc AM~^ is the initial field dependent if 
AM ~ IQ-^Mq, and then B oc AM"^ '^'''Bf is independent of 
tiie initial field after AM ~ 10~^ Mq, until the final stage the field 
B remaining stable at the bottom value of about B ~ 10® G af- 
ter AM ~ 0.2Mo. (2) The bottom field is flie accretion rate (or 
luminosity) dependent as B (x M^^^, which concludes that some 
Z (Atoll) sources with the Eddington luminosity (1% Eddinngton 
luminosity) may possess the magnetic fields of about ~ 3 x 10*^ 
(~ 3 X lO'^) G. Therefore, the existence of the low magnetic field 
(~ 3 X lO'^ G) neutron stars or millisecond pulsars is suggested, 
which needs the confirmation from the future observations. (3) The 
appearance of the bottom field is nothing to do with the initial field 
strength, as shown in Fig.2 and Fig.3. (4) On the NS magneto- 
sphere evolution, Rm oc AM"^ is the initial field dependent if 
AM ~ lO^^Mo , and then Rm oc AM"^ is independent of the 
initial field and the luminosity after AM ~ 10^^ Mq, until the fi- 
nal stage when the magnetosphere radius reaches the star radius. 
(5) The equator magnetic field increases with the accretion and is 
stronger than the initial field strength, which is buried inside the 
crust. Moreover, as known, the occurrence of thermonuclear X-ray 
bursts in LMXBs requires the NSs to have the fields B < 10^" 
G (see Joss & Li 1980), so we may expect the NS in LMXB to 
accrete AM > O.OIM© and to enter into the homogeneous state, 
while the field B forgets the initial field Bo and only remembers 
the bottom field B{. Also, the magnetic field strength distribution 
of flie isolated non-recycled NSs is from ~ 10^° G to ~ lO" G 
(see, e.g., Vranesevic et al 2004; the highest magnetic field can be 
extended to 10^® G if the magnetars are included, see Thompson 
& Duncan 1995), but most MSPs locate in a narrow domain of B 
field distribution, from 10* G to lO" G, which seems to support 
the fact that the NS final magnetic field is not scaled by the initial 
field Bo but scaled by the bottom field Bi. On the characteristic ac- 
creted mass TUB ~ 10~® - 10~^ Mafor the significant deduction 
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Figure 1. The cross section of the accreted neutron star. 9 is the magnetic 
polai' cap angle of neutron star with radius R. H is the thickness of the crust 
and 5H is the expanded immersed depth. Ap and 5Ap are the areas of 
magnetic polar zone and its variation respectively. 
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Figure 2. The magnetic field versus accreted mass diagram with various 
parameter conditions, the Eddington rate L3S and ^36 = O.OIL3S, and 
the initial field strengths from Bg = lO^^ G to Bo = 10^^ G, which are 
indicated in the figure. 
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Figure 3. The magnetosphere radius versus accreted mass diagram with 
various parameter conditions, the Eddington rate L38 and Lag = O.OI-L36, 
and the initial field strengths from Bq = lO^^GtoBo = lO^^G, which 
are indicated in the figure. 



